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Abstract 

In this paper we present a complete classification of the isolated central con- 
figurations of the five-body problem with equal masses. This is accomplished by 
using the polyhedral homotopy method to approximate all the isolated solutions 
of the Albouy-Chenciner equations. The existence of exact solutions, in a neigh- 
borhood of the approximated ones, is then verified using the Krawczyk method. 
Although the Albouy-Chenciner equations for the five-body problem are huge, it 
is possible to solve them in a reasonable amount of time. 



1 Introduction 

The Newtonian «-body problem is the study of the dynamics of « point particles with 
masses m, 6 M + and positions qi E W l (i = 1 , . . . , «), moving according to Newton's 
laws of motion: 

niiinAqj — q ,) 

m M = L 3 !<;<« (D 

¥J r U 

where rjj — \\qi~ qj\\ is the distance between g, and qj. 

In the Newtonian 77-body problem, the simplest possible motions are such that the 
configuration is constant up to rotations and scaling, and each body describes a Keple- 
rian orbit. Only some special configurations of particles are allowed in such motions. 
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Wintner called them central configurations. A configuration (qi,... ,q n ) is called a 
central configuration if and only if there exists aA el such that 

H9j-9G) = l al 3 l<y<« (2) 

¥J K ij 

where qQ = Y,i m i9i/lLi m i i s the center of mass. Equations (0 are invariant under 
rotations, dilatations and translations on the plane. Two central configurations are con- 
sidered equivalent if they are related by these symmetry operations. 

The question of the existence and classification of central configuration is a fas- 
cinating problem that dates back to the 18th century. In 1767, Euler discovered the 
collinear c.c.'s. In 1772 Lagrange proved that, for any three arbitrary masses, the equi- 
lateral triangle is a central configuration. 

For the collinear «-body problem an exact count of the central configurations of « 
bodies was found by Moulton l36l (see also l43l for a modern proof). There is a unique 
collinear relative equilibrium for any ordering of the masses so there are n ! /2 collinear 
equivalence classes. 

The number of planar central configurations of the n-body problem for an arbitrary 
given set of positive masses, has been established only for n = 3: there are always 
five relative equilibria. Two of these are Lagrange's equilateral triangles and the other 
three are collinear c.c. discovered by Euler. Already in the four-body problem there is 
sufficient complexity to prevent a complete classification of the non-collinear relative 
equilibria. In fact, an exact count is known only for the equal masses case (HO an d f° r 
certain cases where some of the masses are assumed sufficiently small l52ll45l . Some 
partial result in the four-body problem with some equal masses are given in l30l[38l l6l. 

Even the finiteness of the central configurations is a very difficult question. This 
conjecture was proposed by Chazy [ lQl and Wintner BP and was listed by Smale as 
problem number 6 on his list of problems for this century 11441 . The finiteness problem 
was settled by Albouy 0] [2 for the case of four equal masses and by Hampton and 
Moeckel lfl6l for the general four body problem. 

Aside from these fundamental results very little else is known in terms of the clas- 
sification of planar c.c.'s for n > 4. Strictly spatial five-body central configurations 
are analogous in some ways to the planar four-body case, and a classification of the 
symmetric ones is given in lfl8l . Planar five-body central configurations present new 
challenges and there are very few results in this area. 

In this paper, we present all the isolated central configurations of the five-body 
problem with equal masses. Such central configurations have been obtained by finding 
all the isolated solutions of the Albouy-Chenciner equations (see Section|2]below, and 
EJUS)) for this problem. Numerical explorations of the problem have been conducted 
by Moeckel fl33l and Ferrario ifTTI . They obtained isolated solutions using root-finding 
routines with a random starting point. However, unfortunately, those numerical ex- 
periments do not guarantee that there are no other solutions, and do not ensure that 
approximate solutions correspond to exact solutions. Our approach is different, and it 
is based on the homotopy continuation method. This is a technique that was proposed 
to find all the isolated solutions of polynomial systems |[T7ll26l . We used the polyhe- 
dral homotopy method to find approximate solutions because it needs to trace fewer 
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curves than other homotopy methods. 

Already for four bodies the Albouy-Chenciner equations are quite large - 12 equa- 
tions in 1 2 variables with total degree 2 , 985 , 984. Moreover, even when the polyhedral 
homotopy method is applied, 82,593 curves need to be traced. Experiments running 
on a 3.2 GHz CPU machine show that the polynomial system solver HOM4PS-2.0 
obtains all the isolated solutions of the four-body problem in 10 minutes, while other 
similar software solvers require more than 10 hours to perform the same task. Consid- 
ering that, in the five-body case, the number of curves that need to be traced reaches 
439,690,761 we use the kernel of HOM4PS-2.0 to find approximated isolated solu- 
tions. We then verify the existence of an exact solution in a small neighborhood of 
the approximate one by using the Matlab toolbox INTLAB, an implementation of the 
interval Krawczyk method |fl9l[35l . 

Note that, while the polyhedral homotopy method finds all the isolated complex 
solutions, this does not exclude the existence of positive-dimensional components of 
the variety defined by the Albouy-Chenciner system in the algebraic torus (C*) 20 = 
(C \ {0}) 20 . Hence, at least in principle, there could still conceivably be undiscovered 
"continua" of positive real solutions. The fact that this is a nontrivial possibility is 
illustrated by Roberts' example of a continuum of real solutions in the Newtonian 5- 
problem with one negative mass |40"1 . 

The paper is organized as follows. The derivation of the Albouy-Chenciner equa- 
tions is shown in Section [2] In Section [3] we introduce the homotopy continuation 
method and several related software packages. The computation results are described 
in Section 21 In Section [5] all the central configurations of the five-body problem with 
equal masses are presented. 

2 Albouy-Chenciner Equations 

The Albouy-Chenciner equations are algebraic equations satisfied by the mutual dis- 
tances nj of every central configuration O [16) . For convenience of the reader, in this 
section, we present a beautiful derivation of the Albouy-Chenciner equations due to 
Hampton and Moeckel lfl6l . 

Multiplying the j-th equation of (0 by mj and summing gives mqa = T!j=i m j<lj> 
where m = Y!}=\ m j- Then, after setting A = mX' , the central configuration equations 
become: 

n 

£ mSij (qi-qj)=0 , l<j<n, (3) 
i=i 

where 

Sij = \+X' (iVj), S tt = 0. (4) 
r U 

Introducing ad xn configuration matrix Q whose columns are the position vectors 
qi we can write equation (f2|i as 

QA = , (5) 
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where A is the n xn matrix with entries: 

A ij = m i S iJ (i^j) Ay = -J\Ay. (6) 

¥J 

The matrix Q can be viewed as representing a linear map from a space of dimension 
n (n being the number of bodies) to the physical space in W 1 where the points are 
located. 

The most important point of the derivation is to replace Q by some quantity which 
is invariant under rotations and translations of the position vectors g, in WL d . Of course 
any such quantity can be expressed as a function of the mutual distances ry. 

Translation of all the positions by a vector u E W 1 transforms Q to Q + uL where 
L is the 1 x n vector whose components are all 1. Thus translation invariance can be 
achieved by restricting the linear map defined by Q to the plane P = {v € R" : Lv = 
v'i + . . . + v„ = 0}. Since the sum of the elements of each column of A is zero it follows 
that A : R" — » P and it restricts to A : P — > P. So QA can be viewed as a linear map of 
PintoR^. 

Rotation invariance is obtained by passing to Gram matrices. For any configuration 
matrix Q, the Gram matrix G = Q'Q is the nxn matrix whose entries are the Euclidean 
inner products qi ■ qj, G is obviously rotation invariant. To maintain translation invari- 
ance, one can view G as representing a symmetric bilinear form j3 (v, w) = v'Gw on P. 
The form /3 on P determines and is determined by the mutual distances ry. To see this, 
note that for any constant ki, adding the vector kjL to the row i of G and the vector k[H 
to the column ; produces a new matrix representing j3. Choosing kj = — \\qi\ 2 shows 
that j3 is represented by the matrix B whose entries are 

Bij =qi-qj-^\qi\ 2 -^\qj\ 2 = —rfj. (7) 

Multiplying both sides of (0 by Q' gives GA = 0. The matrix GA can be viewed as 
representing a (non-symmetric) bilinear form on P, in which case it is permissible to 
replace G by B. Taking the symmetric part gives the Albouy-Chenciner equations for 
central configurations: 

BA+A'B = 0. (8) 

Let e, denote the standard basis vectors in R" and define ejj = e, — ej. Then is 
equivalent to the equations 

e'ijiBA +A t B)e ij = 1 < i < j < n. (9) 

To see this, let y(v, w) = v'Cw be the symmetric bilinear form on P associated to 
the matrix C = BA +A'B. Then (0 means that y(ey, ey) = for 1 < i < j < n. To 
show that y = it suffices to show that y vanishes on the basis en, 2 < i < n of P. By 
the polarization identity 

2y(ei,-,ei;) = y(<?y,ey) - y{e\ h e u ) - y(eiy,eiy), 

this follows from ([9). 
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Equations (0 provide (?) constraints on the (?) mutual distances r, ; - of a central 
configuration. Conversely, it can be shown that if the quantities r, j are mutual distances 
of some configuration in Mr, and if they satisfy (0, then the configuration is central. 
Note that these equations determine the central configurations in all dimensions at once. 

To find the equation explicitly, note that 

Y(eij,eij) = le'ijBAe.j = 2{{BA) n + (BA) n - (M) tj - (BA)jt), 
where (SA)y denotes the entries of the matrix BA. From (O and (|7]i we find 

I m k [Sik{r% - 4 - rfj)+S jk (rl - r% - r?-)] = (10) 

k=l 

for 1 < i < j < n, where Sik and Sjk are given by (@). 

Finally, to eliminate the dilation symmetry, we choose X 1 = —1. 



3 The homotopy continuation method 

Let P(x) =0 be a system of n polynomial equations in n unknowns. Denoting P(x) = 
(pi(x) i ...,p n (x)) and x = (xi, . . . ,x„), we want to find all the isolated solutions of 

!pi (x u ...,x„) =0 
; 
p n (x U ...,X„) =0 

in C". The classical homotopy continuation method 1261 for solving P(x) =0 is to 
find a trivial system Q(x) = (qi(x),... ,q„(x)) and then follow the solution curves in 
the real variable t from t = to t = 1 which make up the solution set of 

H(x,t) = (l-t)cQ(x)+tP(x) = with generic c 6 C\{0}. 

More precisely, all isolated solutions of P(x) = can be found if the system Q(x) = 
0, known as the start system, is chosen properly to satisfy the following three proper- 
ties: 

• Property 0. The solutions of the start system Q(x) =0 are known; 

• Property 1. The solution set of H(x,t) = for < t < 1 consists of a finite 
number of smooth paths, and each of them can be parametrized by t in [0, 1); 

• Property 2. Every isolated solution of H(x, 1) = P(x) = can be reached by 
some path originating at t = 0, that is, the path starts from a solution of the start 
system H(x,0) = Q(x) = 0. 

A typical choice of a start system Q(x) = satisfying Property 0-2 is 

,x n ) = a\x c [ l — b\ 

, Xn ) Q n x n n b n 



CM 



q„(xi 
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Solutions to 
start system 

Q(x) = 



Solutions to 
target system 
P(x) = 




t=0 t=l 

Figure 1 : The solution set of classical homotopy. 



where di,---,d„ are the degrees of polynomials pi(x), . . . ,p„(x) respectively, and 
aj, bj, j = 1,- • • ,n are random complex numbers l24l l25l l26l . The solutions of such 
start system Q(x) = can be explicitly obtained and the total number of solutions is 
d = d\ x di x • • • x d„, which is known as the total degree or the Bezout number of the 
original polynomial system P(x) = 0. We may then find all the isolated solutions of 
P(x) = by following the total degree number of paths originated from solutions of 
the start system Q(x) = 0. But, a great majority of the polynomial systems arising 
in applications have fewer than, and in some cases only a small fraction of, d = d\ x 
• • • x d„ isolated zeros. We call such a system deficient. In this case, many of the 
d\ X • • • x d„ paths will diverge to infinity as t — ► 1 (see Figure [TJ, and those paths 
become extraneous, causing highly wasteful computations. 

In the middle of 1990's, a major computational breakthrough emerged in solving 
deficient polynomial systems. The new method, called the polyhedral homotopy con- 
tinuation method [17|, takes a great advantage of the combinatorial root count, called 
mixed volume, in the Bernshtem's theorem (7), which generally provides a much tighter 
bound for the number of isolated zeros of a polynomial system in the algebraic torus 
(C*)" = (C\{0})". When this method is employed, the number of homotopy paths 
that need to be traced agrees with the mixed volume of the polynomial system. There- 
fore, when the mixed volume of a polynomial system is far less than its total degree, 
the method will greatly reduce the extraneous paths and thereby considerably limit the 
wasteful computations. 

However, this method involves a sometimes costly computation - the mixed cell 
computation. This mixed cell computation can become very costly for large polynomial 
systems. In 2005, a software package, MixedVol fl2l . emerged which led the existing 
codes for the mixed volume computation by a great margin in speed. However, soon 
after MixedVol was published, Mizutani, Takeda, and Kojima l32l developed a more 



6 



efficient algorithm which considerably outperformed MixedVol. Most recently Lee and 
Li ll22l embedded the novel idea of dynamic enumeration of mixed cells of Mizutani, 
Takeda, and Kojima 11321 into MixedVol and a new code, MixedVol-2.0, was produced 
which regained the lead by a substantial margin. 

The software package HOM4PS lfl2l . developed over the years by the group led 
by T. Y. Li at Michigan State University, implemented the polyhedral homotopy con- 
tinuation method for solving polynomial systems. The code was widely considered 
as the most efficient polynomial systems solver. Recently many new curves tracing 
techniques and mixed cell computation algorithms, such as MixedVol-2.0 mentioned 
above, have been included in HOM4PS-2.0. By incorporating these new algorithms 
HOM4PS-2.0 is (as reported in l2T1 ). in general, at least one order of magnitude faster 
than HOM4PS, as well as from one to three orders of magnitude faster than PHCpack 
11471 . PHoM lfl3ll and Bertini 0. Consequently, some large size polynomial systems 
can be solved in a reasonable number of hours. However, as the size of the polyno- 
mial system becomes larger, we need more computing resources to solve the system 
more efficiently. A natural way to attain more computing resources is to parallelize 
the homotopy continuation method which seems to be naturally parallel in the sense 
that each isolated zero is computed independently of the others. The parallel version 
of HOM4PS-2.0, HOM4PS-2.0para, has also been developed. Its excellent scalability 
in numerical results are reported in (27). 

4 Computing central configurations 

The Albouy-Chenciner equations for the five-body problem with equal masses form 
a polynomial system consisting of 20 equations in 20 variables. For finding all the 
isolated solutions in complex field, the polyhedral homotopy method needs to trace 
439,690,761 curves. In order to trace so many curves efficiently, we use the subrou- 
tines in the MPI library (message passing interface 11371 ) to distribute data over multi- 
ple CPUs for parallel computation and we use the kernel of HOM4PS-2.0 to construct 
polyhedral homotopies and trace curves. Employing 32 Itanium2 1.6 GHz CPUs, the 
computation is completed in 140 hours. As a result we obtain 101 , 062, 826 solutions 
in the complex field, 8775 of which are real solutions. Among those real solutions 
only 258 satisfy the physical condition (i.e. ry > for all 1 < i < j < 5). We assume 
that any numerical solution with imaginary parts less than the threshold = 10~ 7 cor- 
responds to a real solution of the Albouy-Chenciner equations. This is reasonable since 
the set of real solutions remains the same when the threshold 6 is chosen between 10~ 13 
and 10~ 3 . In addition, the residuals of the physical solutions are less than 2 x 10~ 14 
and their condition numbers are at most 3.8 x 10 2 , which show that these solutions are 
numerically reliable. 

To guarantee that in a small neighborhood of each numerical solution there is a 
unique exact physical solution, the interval Krawczyk method |[T9l l35l is applied for 
verification. The method is based on the fact: For a smooth function F : W — > W and a 
point x G W, let [x] r Ci* be the interval set centered at x with radius r > 0. Namely, 

[*],- = {yeR":||y-x|L<r}, 



7 



where H-jj^ is the infinity norm. Assume the derivative of F at x, denoted by DF(x), 
is nonsingular, the Krawczyk set of F associated with [x] r is defined as 

K(F,[x] r ) = x-DF(x)- l F{x)+[l-DF(x)- l DF{[x] r )]([x} r -x). 

If the Krawczyk set is contained in the interior of [ x ],., then there exists a unique zero 
of Fin [x],. 

The task of verification is implemented by using the interval arithmetic in INTLAB 
(INTerval LABoratory) PT1 . developed by Siegfried M. Rump at Hamburg University 
of Technology, Germany. In this implementation each numerical solution x is taken as 
the midpoint of the interval set [ x } r with radius r = 10~ 8 . The numerical solutions 
obtained by running HOM4PS-2.0 and the results of verification by using INTLAB are 
available for download at jhttp : 77 hom4ps . math . msu . edu/Nbody . htm . 



5 Central Configurations of the 5 body problem 



In this section we present all the central configurations of the five-body problem with 
equal masses classified according to their dimensions. The dimension of the configura- 
tion was deduced by using Cayley-Menger determinants [4-8 1 . For instance the Cayley- 
Menger determinant for the volume V of the four-dimensional simplex is 
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Configurations of dimension three or less were found by setting the determinant above 
to zero. Configurations of dimension two or less were found by setting to zero the 
determinant above plus the following determinants 
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where, for example, the first of the three determinants above is the Cayley-Menger de- 
terminants for the tetrahedron formed by the masses m\,m2,m?, and 1114. One-dimensional 
configurations were obtained by setting to zero, in addition to the four Cayley-Menger 
determinants above, the ten Cayley-Menger determinants corresponding to all the pos- 
sible triangles formed by a subset of three of the five masses. One of the ten determi- 
nants is 



-16A Z = 
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where A is the area for a plane triangle with side lengths r\2,rys and rn- 

5.1 Collinear and four-dimensional central configurations 

A complete count of the collinear central configurations has been known for almost one 
century. In 1910 Moulton [361] proved that there is a unique central configuration cor- 
responding to each possible ordering of n masses along a line, and thus there are n\/2 
collinear central configurations. In 1970 Smale used tools of the theory of dynamical 
systems to give a modern proof of Moulton's result. 

Our computations show, in agreement with Moulton's result, that there are 60 = 
5!/2 solutions of the Albouy-Chenciner equations. 

The four-dimensional central configurations of five bodies are also well understood. 
Saari l42l proved that the regular n — 1 dimensional simplex is a central configuration of 
n bodies for any value of the masses. In particular the case n = 4 has been well known 
over a century ll23l . The fact that the tetrahedron is the unique spatial central configu- 
ration of four bodies was proved by Pizzetti ||39l . We find one 4-dimensional solution 
of the Albouy-Chenciner equations that, as we expected, is a regular 4-dimensional 
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simplex. 



5.2 Planar Central Configurations 

There are very few results concerning the planar five-body problem. Part of the work 
on the planar five body problem dates back to more than 50 years ago and is due to 
Williams l49l , l50l . More recently, in 2005, Hampton |fT31 found new examples of 
(symmetric) central configurations in the planar five-body problem. Llibre and Mello 
fl28l and Llibre, Mello, and Perez-Chavela |29l found some further central configura- 
tions using the same techniques employed by Hampton. Not much more is known in 
the planar five-body problem 

It has been known for a long time that, in the equal mass case, the regular pentagon 
and the square with one mass at each of its vertices and with the fifth mass at its center 
are central configurations |49l . 

In the case of equal masses we find 147 planar solutions of the Albouy-Chenciner 
equations. Of these, 15 form a square with a particle in the middle (Figure [2] (a)), 12 
form a regular pentagon (Figure 12(b)), 60 form an isosceles triangle with the remain- 
ing particles in the convex hull of the first three (Figure|2](c)), and 60 form an isosceles 
trapezoid with a mass in the convex hull of the first four (Figure|2](d)). The configura- 
tions of Figure|2](c) and Figure|2](d) have been found numerically in 1989 by Moeckel 
ll33l and rediscovered by Ferrario ifTTI . Note that in l49l Williams erroneously states 
that a concave pentagon with three masses at the vertices of a pentagon and the remain- 
ing two in its interior cannot form central configurations. Our results (see figure Figure 
12(c)) show that such statement is incorrect. 

It is interesting to remark that the number of solutions obtained can be explained 
using their symmetry. In the case of the configuration of Figure ©a)) the four masses 
at the vertices of a square and one at the center are subject to an action of the group 
S5 (symmetric group on 5 letters) permuting the masses and hence the set of mutual 
distances r,y. But, the set of mutual distances R for any particular numbering of the 
masses has an isotropy subgroup that is equal to the symmetries of the square (the 
dihedral group D4 of order 8). Hence the set of 55-orbits of configurations ( that is equal 
to the number of solutions of the Albouy-Chenciner equations) of this type should have 
120/8 = 15 elements. Similarly, the pentagon configurations in Figure (|2a)) have 
the dihedral group D5 as isotropy subgroup and there are 120/10 = 12 orbits. The 
configurations in Figure (E](a)) and (E](b)) have C2 (i.e. the cyclic group of order 2) as 
isotropy subgroup and, in this case, there are 120/2 = 60 orbits. 

Another interesting fact is that, as in the four-body problem with equal masses (TJ, 
all the planar central configuration have an axis of symmetry that contains at least one 
mass. 

5.3 Spatial Central Configurations 

In 1908 Brehm J9) found some symmetric spatial central configurations of five bodies. 
Almost one century later Kotsireas and Lazard |fT8l used linear algebra and Grobner 
bases to classify symmetric spatial central configurations of five bodies with equal 
masses. They found four distinct classes of central configurations. Kotsireas and 
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(a) (b) 

Figure 3: Convex spatial central configurations of five bodies with equal masses. 



Lazard conjectured that their list of central configurations was complete. Xia proved, 
for general masses, that there is at least one convex central configuration that is a min- 
imum of the configurational measure l53l . Moeckel used Morse theory, under the 
assumption of nondegenerate critical points, to find 6 more central configurations of 
saddle type |34| . 

We find 50 spatial solutions of the Albouy-Chenciner equations. These correspond 
to four distinct classes of central configurations (two convex and two concave) that 
match the ones of Kotsireas and Lazard, thus proving their conjecture to be correct. 

The first class has 10 solutions. These consist of an equilateral triangle and two 
symmetric points on the axis through the barycenter of the triangle that is orthogonal 
to it (Fig. 0a)). The second one contains 15 solutions. These are pyramidal central 
configurations with a square base (Fig. 0b)). The third class consists of 5 solutions. 
These consist of a regular tetrahedron with a point in its barycenter (Fig. St a)). The 
last class contains 20 solutions. These consist of an equilateral triangle and two points 
lying on the axis through the center of the triangle (Fig. Ub)). 

As in the planar case the number of central configurations can be explained using 
their symmetry. In the case of the configuration of Figure (0 a)) the masses are subject 
to an action of the group S5. But, the set of mutual distances has the group D3 x C2 
(of order 12) as isotropy subgroup . Hence the set of S5 -orbits of configurations of this 
type should have 120/12 = 10 elements and thus there are 10 solutions of the Albouy 
Chenciner equations. The configuration in Figure (0b)) has the dihedral group D4 as 
isotropy subgroup and there are 120/8 = 15 solutions. The configuration in Figure (@] 
(a)) has the tetrahedral group T c i (that has order 24) as isotropy subgroup and hence 
there are 120/24 = 5 solutions. 
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(a) (b) 

Figure 4: Concave spatial central configurations of five bodies with equal masses. 



The configuration in Figure (H](b)) has the dihedral group D3 as isotropy subgroup 
and, in this case, there are 120/6 = 20 solutions. 

According to Moeckel f34l the convex central configurations would fit with the 
Morse theory as follows: 6 saddles plus 9 more =15, 1 minimum plus 9 more = 10. 
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